$1.00 


BASIC Matrix Operations 

















; A” 1963 : 2 I : ‘CATION Or 
W Onn = 1,60 8.20 25, 1( _ | ON OF Ifp 
Mes ae ami] 100, 0 2309 a Families with head ah HEAD, 1968 
] — years and over 
Sa criccrcaccacaesrar eT 
EI oe of School Completed 
Total : ra Some Hi ie 
Families op 2. High gan 
tjo = Schoo! Graduate C nee C 
< ‘O eg 
100.0 DM — 
: : 9 
660 OO 
11.5 es 9.4 Hp 100.0 1 
13.8 17.3 12.6 30 4.0 
a3.2 205 16.8 137 6.5 
25.8 15.0 eo./ 27.4 9.6 
Ale a ve Median ] 2.8 : Ss Z oe 32.7 
CO 1,38 M ncome (Doll 2 “ps 
wicen Income (Dollars. co x, SELECTE? xe CHOY 
or cs : JON: : 
Nite 2. 2,524 Source: Deapy pOPULAT t of total 
deawhite 2 309 Soe ns , 
sidence 
Velen... . e, 1,739 1,613 
ln Central Cities......-...... 676 641 
Outside Central Cities......... 1,064 972 
Nonmetropolitan............... 1,094 1,000 
ther’s Educational Level! 
Less than H. S. Graduate........ 1,153 1,067 
HS. Graduate......... 746 717 
Some College ................-; 306 295 
College Graduate............... 296 279 
mily Income? 
Dader $4000... ............ 501 435 
Sa0eG. (7 o00..:ti‘(‘(‘é‘é#(CN#“(#‘C‘#*(CW 917 855 
7,500-10000................... 521 490 
4 10000-15000........... ....... 508 477 
\ , 15000 and over... ww “ e 169 160 ; 4 
Le ea : 
‘omponents do not add to totals because those who did not report indicated characteristic are ee ee at a A2.\ 58. 
urce: Department of Commerce +380 oo es 
Ll Se ade eet OORT 


O 





BASIC MATRIX OPERATIONS 


The curriculum material and programs in this booklet introduce 
the idea of a matrix. They go on to discuss matrix operations of 
addition, subtraction, multiplication by a scaler, and matrix 
multiplication. The last section covers several contemporary 
applications of matrix manipulations. 


Many versions of BASIC have special "matrix" commands that simplfy 
working with matrices. This booklet does not use these commands. 
There are two reasons for this. First, professional programmers 
do not use them, nor are they available in computer languages 
other than BASIC. Second, you will obtain much more insight into 
the mathematics of matrices by writing programs that manipulate 
them in the same way that the definitions say you should. 


CONTENTS 
Elementary Matrix Operations 3 
Matrix Addition 5 
Matrix Subtraction | 7 
Matrix Multiplication by a Scalar 9 
Storing Matrices in the Computer 11 
Matrix Multiplication L? 
Communication Matrices | Si. 
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& This module introduces the idea of a matrix, and the ele- 


mentary matrix operations of addition (page 5) and subtraction 
(page 7), and multiplication by a scalar (page 9). 


a @ Module #0035 covers matrix multiplication. Module #0045 
discusses several modern applications of matrix multiplication. 


a 
Elementary MATRIX Operations 
The map on the cover of this module shows the streets and 
avenues of a section of the town of Munchin. As can be 
seen, there are employees of Col. Mayer's finger-licking Hot Dog 
Vending Company stationed at strategic intersections in 


the town. Being the business man that he is, Col. Mayer 





is interested in keeping an account of the number of sales 
of each of these vendors. He represents these sales figures 


in an ordered rectangular array as follows: 


SALES FOR MONTH OF APRIL 


qth Avenue 
5 Avenue 
7th Avenue 


12 Avenue 





| This rectangular array is called a matrix (pronounced 
MAY-TRICKS). Each row of this matrix represents locations 
along a certain avenue. Each column of the matrix repre- 
sents locations along a certain street. The entries of the 4 
matrix, which are called elements, represent the number of g 
hot dogs sold by each vendor during the month of April. 

Notice that the elements of the matrix show not only 

the number of sales but also the location of the vendor. 
For instance the entry 877 is the number of hot dogs sold | a 2 


by the vendor located at qth Avenue and 3274 Street. | — 


« @ 


12 Avenue 








1. What is the location of the vendor whose 
sales totaled 4,171 for the month of April? 


ANSWER: | Avenue and Street. 


2. How many hot dogs were sold in April by the 
vendor at 12th avenue and 36th street? 


ANSWER: hot dogs. 


The 0 entries mean no sales because either there is 
no vendor at that corner or the vendor sold 0 hot dogs. 

In the above matrix there are 4 rows (one for each of 
the avenues) and 5 columns (one for each street). The size 
of a matrix is given by the number of rows and columns. In 
this case the matrix is said to have size 4 by 5. [In gen- 


eral, the size of a matrix is given as the number of rows 


by the number of columns. 


MATRIX ADDITION 
Another month has passed and Colonel Mayer has col- 
lected the sales of hot dogs by his vendors for the month 


of May. Here are the numbers: 


_ SALES FOR MONTH OF MAY 










ath Avenue 


5th 


(e [ swe |e 
mf ofan |e [aan 
racer [rise | 0 [sea] 0 


Avenue 


qth Avenue 


th 
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Suppose we now want a table showing the total number of hot 
dogs sold by each vendor for the 2 months of April and May? 
It is probably obvious to you that the way to get this table 
is to add together each vendor's sales for both months. 9 
Doing this for the whole chart is called the operation of | 
matrix addition. 
Before going ahead and doing this, we first will intro- 
duce a more compact notation for matrices,* used by mathe- 
maticians.. We will discard everything except the elements 
of the matrix, enclose them in square brackets, and give 


them a name as follows (capital letters are usually used): 


0 877 0 1,996 
This is the matrix A 
1,652 0 2,008 0 
A= : which describes the 
3,077 2,437 0 2,981 


sales for April 





0 1,002 0 1,890 0 ) 
is j ix M 
1,723 0 2,971 0 1,887 This is the matrix 
M= " 
4,007 2,983 0 3,623 0. which describes the 


sales for May 


3,773 3,076 5, 380 5,179 4,837 





To answer the last question (what are total sales of 
each vendor for April and May) we simply add the matrices 


A and M. This means that each element in A is added to the 


corresponding element in M. - 


*Matrices is the plural of matrix; it is pronounced MAY- 


877 0 1,996 0 
652 0 2,008 0 1,765 
077 2,437 0 2,981 0 


1,002 0 1,890 0 
723 0 2,971 0 1,887 
007 2,983 0 3,623 0 





0+0 877+1,002 0+0 1,996+1,890 0+0 
The matrix ,6524+1,723 0+0 2,008+2,971 0+0 1, 765+1,887 
A+M 3,077+4,007 2,437+2,983 0+0 2,981+3,623 0+0 
3,618+3,773 3,051+3,076 4,1714+5,380 4,632+5,179 3,987+4,837 
0 1,879 0 3,886 0 
Spada 0 4,979 0 3,652 
| KKKKKKKKKE 
. *EXERCISE* 
. KEREKEREKK 
z 
Fill in the missing entries in the sum matrix shown 
above. (Hard work? Cheer-up. You'll learn to use the com- 
puter for this purpose soon.) 
MATRIX SUBTRACTION 
f Colonel Mayer wants to know the difference in 
a 


hot dog sales between May and April for each vendor so that 
he can tell which location is showing the greatest boom 


(increase) in business. 


The answer to this problem involves the matrix opera- 
tion of subtraction. Matrix subtraction is similar to 
matrix addition, except that each element of A will be 
subtracted from the corresponding element of M. Thus, | 
> 
0 1,002 0 1,890 0 
1,723 0 2,971 0 1,887 
M-A= 
4,007 2,983 0 3,623 0 
Splis 3,076 5, 380 5,179 4,837 
0 877 0 1,996 0 
1,652 0 2,008 0 1,765 
= ; Sa 
3,077 2,437 0 2,981 0 
3,618 2,051 4,171 4,632 3,087 
1,002- 877 0-0 1,890-1,996 0-0 
,723-1,653 0-0 2,971-2,008 0-0 1,887-1,765 
M-A = 
4,007-3,077 2,983-2,437 0-0 3,623-2,981 0-0 
| % 
0 125 0 -106 0 , | 
70 0 963 0 122) .- 
930 0 0 
155 850 
XY 
KkkeKkkKkkkkRKK Tia 
*EXERCISE* 
KEEREREKRES 


1. Complete the above "difference" matrix. 


In performing the matrix operations of addition and 
subtraction, notice that the 2 matrices being added or sub- 


tracted must be of the same size. 


MATRIX MULTIPLICATION BY A SCALAR 

With the approach of the summer months and the baseball 
season, Col. Mayer is Srediovina a 20% increase over the 
May sales for the month of June. What number of hot dogs 
does he predict each vendor to sell during the month of 
June? 

The answer to this question involves the matrix oper- 
ation of multiplication by a scalar. This means that each 
element in a matrix will be multiplied by the same number 
which is called a scalar. The scalar multiplier for our 
example is 1.20. We can say that each element in June is 
predicted to be 1.2 times as big as the corresponding ele- 


ment in May. The compact way of saying this with matrices is 











1.200 1.20*1,002 1.200 1.20* 1,890 1.20x0 
1.20*1,723 1.20*0 
J= 
1.20*4,007 1.20*0 1.20*0 
1.20*2,773 1.20*3,076 1.20«5,380 1.20*5,179 1.20*4,837 
0 0 0 
ED A 0 eS 
l= 
0 0 


- 10 - 


* 
x (What, again!) 
* 


Complete the matrix J=1.20*M shown on the previous page. 
(If you get tired, read on and find out how to use the com- 


puter to do this.) 


Using the Computer to Add and Subtract Matrices and to Multipl 
Matrices by a Scalar 


n IMPORTANT NOTE? n Many versions of BASIC have special 
"matrix" commands that simplify working with matrices. 
However in this module we will not use these commands. 
There are two reasons for this. First, professional 
programmers do not use them, nor are they available in 
computer languages other than BASIC. Second, you will 
obtain much more insight into the mathematics of matrices 


by writing programs that manipulate them in the same 


way that the definitions say you should. 


This module assumes that you understand what a subscripted 


variable is, and how such variables are used in BASIC. 







READ ABOUT SUB- 
SCRIPTED VARIABLES 
AND THE DIMENSION 
STATEMENT IN 
"“Populution: A 
SELF-TEACHING 
BASIC PRIMER" 





S THIS A COR- 
RECT ASSUMPTION 













CARRY ON! 
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i STORING A MATRIX IN THE COMPUTER 
Since an element of a matrix is located by its row 

number and column number (see page 2), the best way to store 

matrix elements is to use variables (locations) which have 

= two subscripts. Then we can store the element which is 

o in row i and column j in the location called A(I,J). For 
example, A(1,2) will be the variable name we give to the 
location where we store the matrix element in the first row, 


second column of A. 


INPUT OF MATRIX ELEMENTS 
The program segment shown below allows you to input the 
elements of a 2 by 3 (2x3) matrix. It expects you to input 


the numbers left to right, row by row. 






ene 
105 DIM A(2,3) After you run this program, the 
110 FOR I=1 TO 2 memory of the computer looks. 
115 FOR J=1 TO 3 like the following: 
120 INPUT A(I,J) 
125 NEXT J 
130 NEXT I ) 
1000 END A(1,1) A(1,2) A(1,3) 
RUN 
212 12 “2. 14 
2-2 
214 A(2,1) A (2,2) A (2,3) 

' ?0 

© 233 0 33 19 

?19 





PRINTING OUT ELEMENTS 


To print out the above matrix elements row by row we can 
add the following statements to the above program: 


910 FOR I=1 TO 2 


| 915 FOR J=1 TO 3 The semi-colon forces all the elements 
920 PRINT A(I,J); controlled by the inner FOR Loop to be 
‘ 925 NEXT J printed on one line. 


a 930 PRINT (cr) The extra PRINT statement forces the com- 
935 NEXT I puter to print the next row on a new line. 


a A as 


Program (on p. 4). 


Matrix Addition: C=A+B 


1 DIM A(4,5), B(4,5), C(4,5) S 

5 PR. "THIS PROGRAM ADDS TWO 4 by 5 MATRICES" 
10 PR. "INPUT THE ELEMENTS OF MATRIX A" 

15 FOR I=1 TO 4 


20 FOR J=1 TO 5 
25 INPUT A(I,J) 


30 NEXT J 
35 NEXT I 


40 PR. "NOW INPUT THE ELEMENTS OF MATRIX B" 
45 FOR I=l1 TO 4 


50 FOR J=1 TO 5 
55 INPUT B(I,J) 


Input Matrix B 


60 NEXT J. 

65 NEXT I 

70 FOR I=1 TO 4 This segment of the program 

75 FOR J=1 TO 5 adds each element of B to the 
80 C(I,J)= A(I,J)+B(I,J) corresponding element of A, 

85 NEXT J and puts the sum in the cor- 

90 NEXT responding position of Matrix C 


95 PR. "THE SUM MATRIX IS: 
100 FOR I=1 TO 4 
105 FOR J=1 TO 5 


110 PRINT C(I,J); This prints out the Matrix C 
115 NEXT J which contains the sum elements 
120 PRINT 

125 NEXT I 

130 END 


Write a program to calculate C=A+B for matrices that 
can have any size up to and including 10 rows and 10 columns. 


HINTS: 
a. 


b. 


Dimension A, B, and C as 10 by 10 


Allow the user to input the number of ROWS = M 


Allow the user to input the number of COLUMNS 


Change all the FOR loops to go TO M or TO N 


Input Matrix A 
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i A 


CHALLENGE PROGRAM #2 
PPP Poo PoP P PP PDD??? POI? 


Test the above program on your own 10 x 10 matrices. 


ALLENGE PROGRAM #3 
Ddd3 2222722222222 P 2??? 


VG 


fs 


Write a program which inputs a Matrix A and inputs a 
scalar factor J. The program should output the Matrix D 


which is the scalar product of J times M, i.e. 


MATRIX SCALAR — MATRIX 


THIS IS A 
SINGLE NUMBER 


THESE ARE ARRAYS 
OF NUMBERS 











Test this program on the Hot Dog Problem on page 6. 


¢ 


‘es 
CHALLENGE PROGRAM #4 
WUE EN EEEUAEE EE SEE EE 


(Confidential Note: We will throw in some technical 
language in this problem which will make it appear compli- 
cated. But the mathematics (and computing) is no harder 


than that used in the Hot Dog problem.) 


= JA om 


The chart below shows the observed frequency values 
found in a series of pathologically verified cases of ulcer 


sease * 
and cancer of the stomach: 


aes Achlorhydria | Hypochlorhydria Hyperchlorhydria 
CS A 


The expected frequency values are given in the table E 







below: 


BP Achlorhydria | Hypochlorhydria Hyperchlorhydria 


These tables can also be represented as matrices: 










= 3 i 35 9 | (observed values) 
22 2 6 0 
E= 16 a2 20% Sait (expected values) 
9 3.2 14.6 3.2 


For statistical purposes we would like a table showing the 
differences between the observed and expected values. This 
means that you should modify your program to output the 
matrix F-E. HINT: As you might guess, all you have to do 
is modify the matrix addition program on page 9 so that line 
80 uses a 'minus sign' (-). (Also change line 95.) 


—— 
See the book "Facts from Figures" by M. J. Horoney, pp. 246-270. 


og. ae 





Modify the above program so that a matrix of relative 
differences between observed and expected values is printed 


out, followed by a matrix of percentage differences. 


Example for one element: (let's take the first row, second 
column) 


¥ (1,2) 


7 (observed) 


E(1,2) 5.8 (expected) 


The relative difference by which the observed value differs 
from the expected value for these particular elements is 


calculated by dividing the difference between F and E by E: 


(7-5.8)/5.8 = 0.20689 = RELATIVE DIFFERENCE 


= The percentage difference is obtained by multiplying 


the relative difference by 100: 
0.20689%*100 = 20.689% = PERCENTAGE DIFFERENCE 


Now go ahead and do this for all elements in the matrices, 


using the computer as suggested on the next page. 


1 2a = 


HINT 1: Modify line 80 (p. 9) so that the Matrix C 


will contain relative differences. 
C(I,J) = (F(I,J)-E(I,J))/E(I,J) 


HINT 2: For percentage differences do a scalar multi- 
plication by 100. You can use the Matrix C "over again" 
(after printing the relative difference) by letting C(I,J) = 


100*C(I,J) within a FOR loop. 


HINT:3: Since you will want to print out two matrices 
(relative differences and percentage differences), but since 
they both will be called C (see HINT 2 above), it would be 
a good idea to put the print statements in a sub-routine as 


follows: 


INPUT E and F 


FILL THE MATRIX C 
WITH RELATIVE DIFFERENCES (See HINT 1) 
GOSUB 800 


FILL THE MATRIX C 
WITH PERCENTAGE DIFFERENCES (See HINT 2) 
GOSUB 800 


SUBROUTINE TO PRINT C 


RETURN 





800 







ae ome 


— Matrix Multiplication 
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PROBLEM 
x . 
A B | means Matrvtvx 
. Multtpltcation!!? 
5 4 2 27 2 9 
39 81 * 15 3/=]? ? 
621 6 4 ? ? 
that means I have to 
SOLUTION use real number multt- 
pltecatton 
(+) and 
addttton 
(+) tn 


just the 
rtght way. 
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PROJECT SOLO / Dept. of Computer Science / Univ. of Pxttsburgh(15213) 
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MATRIX MULTIPLICATION 


INTRODUCTION 
The need to develop a compact notation for manipulating 


systems of linear equations led 19th century mathematicians 





to the discovery and use of matrices. (Deep question: Ane 
mathematical concepts discovered or invented?) For the def- 
inition of a matrix and a discussion of matrix addition and 
multiplication of matrices by scalars see the module "ELE- 
MENTRAY MATRIX OPERATIONS". The present module is a presen- 
tation of another important matrix operation--matrix 
multiplication. 

Multiplication of matrices was first suggested by re- 
search in the theory of linear equations, but it has turned 
out to be an important concept for many other applications, 
including mutivariate statistics (used by all soctak sceken- 
41sts), analysis of communications within a group (see 
module 0045 "COMMUNICATION MATRICES"), electrical network 


problems, and mathematical models in psychology and economics. 


Let's start with a SPECIAL CASE: Multiplying A Matrarx 
by a Vector. | & 


Consider the following set of three linear equations 


in three unknowns: 


4x, + 1X5 + 12x, = 95 
3X} Se 2X5 + 10x3 = 58 
+ 8x + 5x = 121 

Ax, 3 3 _ 


Se eet 





GZ 


= 7S. 


Matrix multiplication is defined in such a way that the 


above set of equations can be expressed in compact matrix 


form as follows: 


Where: 


NOTE: 


A *xXs B 
1. "*" is the symbol for matrix multiplication 


2. A is the matrix of coefficients 


4 7 12 
A = 3 2 10 
11 8 5) 


3. B is the column vector of constant terms. 
A cokumn vector 48 a matrrcx wath only one 


column. 

95 

B = 58 
121 

4. X is the column vector of unknowns. 

a] 

X = Xo 
3 


Sinee X 45 a matrix 4t45 eLements should have two 
Subsenripts--the first for the now and the second 
for the cokLumn--that is, strictly speaking we 
Should write: 


*31 
But sinee there 4b onky one cokumn there should be 
no confusion 4§ we drop the second (column) subscript. 


=> Bh i= 


That is, matrix multiplication is defined so that: 


4 7 12 X41 95 
3 2£ AO . Xo = 58 
ll 8 6. X3 121 


is a true statement. 


But notice that 


95 = 4x, + 1X5 + 12x, 
58 = 3x4 + 2x5 + 10x, 
121 = 11x, + 8x, + 9X3 


hence, what we have said is that 


4 7 12 xy 4x, + 7X + 12x, 
3 2 10 * Xo = 3X} + 2X. 4 10x, 
ll 8 5 | X3 11x, + 8x. + 9X3 


is a true statement by the definition of matrix multiplication. 


In general, if 


A = aoj ao ao3 3 AOWS 
497 Sap “a9 
then Ail ai a3 Xy 
* = 
A*X Ao, 222 493 re Xo 3 AOWS 

45, 892 #9 *3 

1 column 

471%, FT 41 o%2 * 343%3 


a51%*] -+ Ao9X5 + 253% 3 NOWS 
25174 @g5%o F 2955 


NOTE: This has onky 1 column 
one column. 


C ) 


= ot im 


QU Notice that the product of a 3 x 3 matrix and a 3 x 1 


matrrx (cokumn vector) ib a 3 x 1 matnix (column vector) 


which is formed as follows: 


Multiply each number in the {415 now of A by the cor- 





responding number ({i145t by first, Second by second, third 





by third) in the §4445t column of X [in this case X only 
has one column) and add these products. This gives the 


number in the first row and first column of A * x, 


now 1 column | 
how 1 


492 9172 913] * = [*]} = | 212%1 * 212*2 * 213%*3 
a. 
e @ ® an. 6 e e 
column 1 
Do the same for the second row of A and the {inrst, and 


only, column of X to determine the number in the second 


row and first column of A *X. That is: 


1 how 2 column 1 
how 2 aoj ao a5 3 * Xo = 251%] + AooX5 + 253%3 
e e e x 3 , e 
cofumn 1] 


Finally the number in the thind row and the {inst column 





of A * X is computed as follows: 


2 ® ° 9 aay e 
e e e Xo e 
wow S$ | 3, 432 333] * |¥3] = |232%1 * 232%2 + 295%3 


column 7 now 3 column 1 


= 39 = 


To repeat the final result: 


ee, ti ee 
451 912 333 ay 9 %4t Ay oXat 4] 4X3 
$91 “95 . “55 * |Xo) = 291%4t AgoX2t a53%; 
231° 839 “S95 ac 437% + AzoXot a3 5%;3 





NOTE: To keep track of what you ane doing, it is sometimes 


~ helpful to point to the elements that are being . 
multiplied. 
For example, when you are computing the number in row 2 
column 1 of A * X above, 
Your Lest hand points to: Your rAtght hand points to: 
“93, “99 “og 
Sie 















221%1" a 


eR HH Li IP aR EK AB 


22%2* 22 ya)4 


etc, 


Try it--you'll like it! 


— ¥e = 


a EXERCISE 1A A* xX =B 


A program in BASIC which multiplies 3 x 3 matrix* by 


a 3 x 1 column vector might be organized as follows: 





a 
SET ELEMENTS OF B = See module 107 on 
READ ELEMENTS OF A "ELementary Matarix 
AND X Operations" 





FOR I iL TO 3 
FOR J 1 TO 3 
B(I) = B(I) + A(I,J) * X(J) 


NEXT J 
NEXT I 


— PRINT 
B 





Write the complete program and try it with the 


following matrices: 


2 | 3 1 7 -] 1 0O 0 20 
N = 3, 9 5 2 * 6 N = 3, 0 1 0 * | 20 
2 4 8 4 0 O 1 20 
Also multiply these by hand as well as by using your program. 
‘eet 


* 3 x 3 is "pronounced" 3 by 3; the x does not mean multiplication. 


as 94 & 


EXERCISE 1B 


In the above discussion there was nothing special about 


the number 3. The more general definition for multiplying a 


N x N matrix by a N x 1 matrix (column vector) 


44, F412 °** Fin ay 457%] TF 4, 9%*9 
491 992 °° Bon] « | *2) = 1221%1 * 222*%2 
an an2 a ann xn 11 ¥ an2*%2 


-- « 


+ a 
+ a 
+ a 


Write a BASIC program which reads a number N, aN x N 


matrix, and a N x l matrix, computes the matrix product, and 


prints the result. 


Programming Note: A feature of BASIC which 4&8 useful an 


dealing with matrices is SUBCRIPTED VARIABLES. 


See YOur 


BASIC reference manual. You may also want to review 


FOR-NEXT Loops. You wikhk have to decide on a maxdmum 


possible vakue for N and write a DIM statement accordingly. 


EXERCISE 1C 


Test your program on the following data: 
(a) (b) 


1 2 3 4 5 1 21 73 
6 7 8 9 10 2 18 98 
11 12 13 14 #215 n 3 26 47 
16 17 #18 19 20 4 22 27 


21 22 23 24 25 


Ul 


84 
41 
93 
68 


cc 


= 95 = 


The Genenak Case: Mone Than One Column 


Suppose in the product A * X = B, A is a 3 x 3 matrix 
and X is a 3 x 2 matrix. Matrix multiplication is defined 
so that B will be a 3 x 2 matrix defined such that: 


Each column of the product matrix B is the 
column vector that you get by multiplying 
A by the corresponding column of X. 


This is illustrated very well on the cover of this module 
and it would be worth your time to study it carefully. 


For those who like symbols, subscript, etc. 


Si, “Za "13 “1 “12 
5 ROWS = Van, Ana 823] * 1X2) X22] = 
“ai. “go “gg “31 “92 
3 columns 2 columns 
Sy" A494 © Aas, | Shae A 1e%oo * ByaX a2 
ee Goa T Sod" oq * Bga%ay Bape F Syg%oa * So9%39 
a3 “Y Soa%s1 © Agg%s1 “= gitMaa  *39%2o * Aa9% a2 


2 columns 


Matrix multiplication is defined for matrices A and X 
whenever the number of columns of A is equal to the number 
of rows of X. Another way to say this is that any N x K 
matrix can be multiplied by any K x M matrix. The product 


matrix will be N x M matrix. 


The rule to remember is: [N x K] * [K x M] = [N x M] 


must be the same 


sit Gs 


Matrices which can be multiplied are called conformable. _# 
Notice that it may be possible to multiply A * X but not 


possible to multiply xX *« A. (Can you think of an example?) 


EXERCISE 2A 





Write a program which inputs three numbers N, K, M, a 
matrix A with N rows and K columns, and another matrix B with & 
K rows and M columns. The program should output the matrix 
product A * B. 

EXERCISE 2 B 
Try your program on those products below which are conformable. 


a) N= 3, K=3, M=2 
5 4 2 27 


3 9 8] * 3 35 


ll 
ae) 


621 6 4 Seo” 


b) 


NO 
© il 
Nm Ww 

= 
A 
| 
ke WwW 
ay 
oO 
= 
(>) 
ll 
WwW 


~] 
10.0) 
© 
© 


EXTRA: An element I (in a Set with a multiplication operation 


defined on it) such that I * A =A * IT =A for any A 4b 


called a multiplicative identity. What is the multiplicative @ 
Kdentifty for Nx N (Square) matrices? What ts the additive 


Adentrty? 
Answers for N = 2 





MULTIPLICATIVE ADDITIVE 
Identity 0 Identzty 


 ) 





a DF. mw 


eee |S -3 1 0 2 4 =2 
-“l 3-1 3 6 4-5 12 8 3 
* = 2 
2 0 5 6 7 6-1-3 #=5 9 
“7 1 2° 3 Oo -1l1 5.°9 -7 


The above two matrices are not conformable for multi- 
plication in the reverse oder. Hence, multiplication is not 
commutative for every pair of matrices. (That is, it is not 
true that A * B=B* A fon ahh matrices A and B.) But 
what if A and B are N x N (Square) matrices? Do you think 


that matrix multiplication is commutative for square matrices? 


d) N=2, K= 2, M=2 


What can you conclude 
from (d) and (e)? 


ss 
No 
KH Oo 
| ee 
+ 
bP 
NO J 
es 
il 
J 


£) N= K=M= 3 
1.000000 0.500000 0.333333 5 -36 30 
0.500000 0.333333 0.250000 * |-36 192 -180 
0.333333 0.250000 0.200000 30 -180 sail 
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LISTING OF A GENERAL PROGRAM TO MULTIPLY MATRICES 


SSaAssssss 


PRINT “THIS IS A PROGRAM TO MULTIPLY 2 MATRICES.” 
PRINT “THE MATRICES MUST BE CONFORMABLE & OF” 
PRINT “SIZE 10 BY 10 OR LESS-” 

PRINT “INPUT THE DIMENSIONS OF THE MATRIX TO BE MULTIPLIEDe”™ 
PRINT “WHAT IS THE NUMBER OF ROWS?” 

INPUT RI 

PRINT “WHAT IS THE NUMBER OF COLUMNS?” 

INPUT Cl 

PRINT “INPUT THE DIMENSIONS OF THE MULTIPLIER MATRIXe™ 


PRINT “WHAT IS THE NUMBER OF ROWS?” 

INPUT R2 

IF R2@2Ci THEN 150 

PRINT “THE MATRICES ARE NOT CONFORMABLE.“ 
GOTO 90 

PRINT “WHAT IS THE NUMBER OF COLUMNS?” 
INPUT C2 

PRINT “INPUT THE MATRIX TO BE MULTIPLED:” 
PRINT “INPUT BY ROWSs ONE ELEMENT AT A TIME” 
FOR I=1 TO Ri 

FOR Js! TO Cl 

INPUT ACI »sJ) 

NEXT J 

NEXT I 

PRINT “INPUT THE MATRIX THAT WILL BE THE MULTIPLIER.” 
PRINT “INPUT BY ROWSs ONE ELEMENT AT A TIME” 
FGR I2#1 TO R2 

FOR J2=1 TO Ce 

INPUT BCI.sJ) 

NEXT J 

NEXT I 

FOR I=#1 TO Ril 

FOR J=1 TO C2 

LET CC€1s2J)=0 

NEXT J 

NEXT I 

FOR Kz1 TO C2 

FOR I2#1 TO Ri 

FOR J=1 TO R2 

LET CC1seK) #CC Is KI +ACIs J) *BC Js K) 

NEXT J 

NEXT I 

NEXT K 

PRINT “THE PRODUCT MATRIX IS®* 

PRINT 


4A0FOR I#1 TU Ri 


FOR Jel TO C2 

PRINT CCIsJ)3 

NEXT J 

PRINT 

NEXT I 

PRINT 

PRINT “DO YOU WISH TO DO ANOTHER MULTIPLICATION?” 
PRINT “IF SO TYPE A 1 OTHERWISE TYPE A 2e"” 
INPUT Y 

IF Y=! THEN 40 

END 
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SAMPLE RUN OF /MAT/ USING THE DATA FROM THE COVER PICTURE 


>RUN 

THIS IS A PROGRAM TO MULTIPLY 2 MATRICES.~ 

THE MATRICES MUST BE COMPATIBLE AND OF 

SIZE 10 BY 10 OR LESSe 

INPUT THE DIMENSIONS OF THE MATRIX TO BE MULTIPLIED.- 
WHAT IS THE NUMBER OF ROWS? 

2 3 

WHAT IS THE NUMBER OF COLUMNS? 

? 3 

INPUT THE DIMENSIONS OF THE MULTIPLIER MATRIXe 
WHAT 3S THE NUMBER OF ROWS? 

2? 3 

WHAT IS THE NUMBER OF COLUMNS? 

22 

INPUT THE MATRIX TO BE MULTIPLEDe 

INPUT BY ROWSs ONE ELEMENT AT A TIMEs 


? 5 

t 

? 2 

? 3 

29 

? 8 

? 6 

22 

? 3 

INPUT THE MATRIX THAT WILL BE THE MULTIPLIER. 
INPUT BY ROWS» ONE ELEMENT AT A TIME. 
? 2 

? 7 

25 

? 3 

? 6 

? 4 | 

THE PRODUCT MATRIX IS 
42 55 

99 80 

28 Ss2 


DO YOU WISH TO DO ANOTHER MULTIPLICATION? 
IF SO TYPE A 1 OTHERWISE TYPE A 2e 
2? 2 


PEXIT 


~ 30 = 





a 
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— Communication Matrices 


| PROJECT SOLO 
Depantment of Computer Science, University of Prttsburgh 
Pattsburgh, Pennsylvania (15213) Module 0045 


(FEATURING THE CASE OF THE MUTATED VIRUS ON PLANET X!) 





@ This module shows the application of matrices in studying 
i communication patterns. Problems related to airline ticketing, 
message delivery, and communicable diseases will be presented. 

& The module "Matrix Multiplication" is a prerequisite for do- 


ing the problems in this module. 





The Airline Ticketing Problem: 


ay ee 


EAGLE Airlines, a little known but ambitious "scheduled" 


commuter service, lists the following direct flights: 


Eichelbergertown 
Eichelbergertown 
Eichelbergertown 
Applebachsville 
Applebachsville 
Lackawaxen 
Lackawaxen 
Lackawaxen 
Satchelville 
Satchelville 
Salladasburg 
Klingerstown 


Klingerstown 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


to 


Salladasburg 
Klingerstown 
Lackawaxen 

Dead Man's Gulch 
Klingerstown 
Applebachsville 
Klingerstown 
Salladasburg 
Klingerstown 
Eichelbergertown 
Eichelbergertown 
Satchelville 


Salladasburg 


An interesting question we can ask is the following: 


For any given starting city, how many different 
ways can you fly to each of the cities (including 
the starting city) via EAGLE Airlines? 


(We are not concerned with how long a passenger may have to wait 


for a connection in this problem). 


a AS 


Diagrammatically we can represent EAGLE's flight schedule 


as shown below: 


@DEAD MAN'S GULCH 





APPLEBACHSVILLE 
- 
LACKAWAXEN 
KLINGERSTOWN 
SALLADAS BURG 
SATCHELVILLE EICHELBERGERTSON 
A diagram such as the above is called a directed graph. The 
direction of the arrows represents a flight from one town to an- 
other. 
(NOTE: We're assuming that there are no flights scheduled that take 
off, circle overhead, then land in the very same town! Such a con- 
% nection would be represented as: 
2 Such connections are useful in 


some applications. After finishing KLINGERSTOWN 
this module, see if you can think of 
any such applications. Also check the modules on "Finite State 


Automata." ) 


Another way to represent the EAGLE flight routes is with 
a rectangular array: 


. TO: EICHELBER-| APPLEBACHS- | LACKA- | SATCHEL- | SALLADAS-+ KLINGERS-| DEAD MAN'S 
FROM GERTOWN VILLE WAXEN VILLE BURG TOWN GULCH 


EICHELBER- 
pena NO NO YES NO YES YES NO 


APPLEBACHS | 
VILLE Ne NO | NO NO NO | YES | YES 


LACKA- 
WAXEN 


NO | 


SATCHEL- 
VILLE NO | 


SALLADAS- : | 
BURG NO : 


KLINGERS- 


TOWN YES NO 





DEAD MAN'S 
GULCH 


Either "YES" or "NO" is put in each box in answer to the ques- 


tion "Is there a direct flight from (row city) to (column city)?" 


A part of the table has been filled in. Complete the table 

by looking at EAGLES's schedule on page 2, placing a YES or NO in each space. 
When using the computer to help solve Problem 1, the table 

above is represented as a matrix. The rows and columns of the 

matrix correspond to the rows and columns of the table. The 

elements of the matrix are 0 to represent NO and 1 to represent 


YES. 





Ft s 








=, BE. e 


The matrix representation would be: 


Soo OC OrF.O: 
ooo 0 CO Oo 
SoFOOO oO oO 
SOorFOOFrFO RF 
cCOOrF FF 


By representing the flight schedule as this binary matrix 
we can extract further information about the schedule. The 
matrix operations of addition and multiplication are used in 


the investigation. 


In the binary matrix the non-zero entries represent a direct 
connection from one town to another. This direct connection will 
be called a l-path, or a path of length-1l. By squaring this matrix 
(raisingit to the 2nd power i.e. calculating A*A) we produce a 
matrix whose entries are the number of 2-paths, or paths of length 
2, from one town to another. Such a path of length 2 requires 2 


direct connection flights to get to the destination. 


Squaring the matrix A gives: 


4) 
oO 
rh 
Q 
()) 
Qu 
iQ 


bike. bawO. + 2,2 . 0) (NOTE: A*A was calculated 
aN using the method explained 
0 0 0 : | 1 2 \ in the module "MATRIX 
1 0 0 2 2 2 1} \ MULTIPLICATION". ) 
a?caean = |(0° O° Bb (1 °-2: 1 % 
0 Oo 1 O +t 1 0 First row of A?” 
2 0 0 0 0 1 0 
0 0 0 0 0 0 0 


— 


WHAT DOES THIS NEW MATRIX TELL US? = 


Look at the first row. The numbers there represent the 
flight originating from Eichelbergertown. There is one 2-path to 
each of the towns: Eichelbergertown, Applebachsville, Satchelville oe 
and Klingerstown, as follows: 

q 

a. Eichelbergertown to Salladsburg to Eichelbergertown 
b. Eichelbergertown to Lackawaxen to Applebachsville 
GS. Eichelbergertown to Klingerstown to Satchelville 


d. Eichelbergertown to lLackawaxen to Klingerstown 


There are two 2-paths to Salladsburg: 

Eichelbergertown to lLackawaxen to Salladsburg 

Eichelbergertown to Klingerstown to Salladsburg 
Ney 


There are no 2-paths to Lackawaxen (f) or to Dead Man's Gulch (gy 


What do you think the entries of the matrix which results from 


cubing the original matrix represent? 


We willuse the mathematical notation A* to mean A*A, and A? 


to mean A*A*A; Note that A? also = A7*A, 


In our example with seven cities (in general n cities) it + 
should be “obvious” that a passenger can eet from one city to 
another city with six (in general n-1) or less direct flights? 
If we also include round trips (that is, trips which return the 
passenger to his starting city), we should count the paths of 
length nor less. Thus in Problem 1 we're interested in counting rn 
all the paths of length n or less that connect any given starting -_. 


City with any other city. 


* Assuming he can get there at all. 


a en 


One method of obtaining this information is to look at the 
diagram on page 3 and count all of the different ways. However, 
it is very easy to miss some paths. A more sophisticated method 


is to sum all of the power matrices (lst power, 2nd power, etc.) 





from the lst power to the nth power, where n is the number 
e of cities (and therefore the number of rows (or columns) in the 
binary matrix A). [If you don't remember how to add matrices, work 
through the module "ELEMENTARY MATRIX OPERATIONS". ] 

In the mathematical notation this is 


l 
% «o 


S=A!+a? +ai+... + ae” 
where S is obtained by adding the matrices on the right. The 
entries in the summation matrix are the total number of paths of 


any length (1 ton) connecting any two towns. 


PROBLEMS FOR COMPUTER SOLUTION 


Write a computer program which finds the number of different 


Problem 1 


ways in which a person can fly from a givenstarting city to each 


of the other cities (including the starting city) via EAGLE Airlines. 


Sample Solution 


: On page 8 we will show a flow chart describing one aeturien of 
this problem; on page 9 a listing of a BASIC program that implements 
this flow chart is given. This program doen not use the matrix 
functions available in most versions of BASIC in order to show you 


how a professional programmer manipulates matrices. 


ee Se 





Flow Chart, Problem l Se 
- AY | 
40-65 2 yh | 
a 2 
70 
er er on ee ee 7 
75-95 B* Ak | : 1s a cc aie GEA. | 
—Teaa  # a = x”. | 
Lb — a 
| Adds the Nth power matrix | 
100-120 o + CO eee | to the sum matrix. 1 
gen 4 A° @ ca, ww, | 
Rial SCeE. A Wrecc mo Nines OE REI ATTEN - 
125 N+ 1 esa” 
165-195 


130 IF N> 7 GoTo 165 | Le J 165 PRINT MAT S 
~™ 
~ 





| 
B now becomes A | 
for the value of N l 
calculated in step 6. | 
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Sample Program 


5 DIM A(7,7), B(7,7), C(7,7), S(7,7) 

10 PR. "INPUT THE CONNECTION MATRIX ROW BY ROW" 
NY 15 FOR I= 1 TO 7 

20 FOR J= 170 7 

25 INPUT A(I,J) 






30 NEXT J , 
35 NEXT I The line numbers used here 
correspond to those shown on 
40 FOR I= 1 TO 7 the flow chart of page 8. 
45 FOR J= 1 TO 7 
50 S(I,J) = A(I,J) 
; 55 B(I,J) = A(I,J) 
v 60 NEXT J 
65 NEXT I 
70 LET N= 2 
75 FOR I= 1 TO 7 
80 FOR J= 1 TO 7 
*85 C(I,J) = B(I,1)*A(1,J)+B(I,2) *A(2,J)+B(I,3)*A(3,J3)+B(L,4) *A(4,J) 
+B{I,5)*A(5,J)+B(I,6) *A(6,J)+B(I,7) *A(7,d) 
90 NEXT J 
95 NEXT I 
100 FOR I= 1 TO 7 
105 FOR J= 1 TO 7 
110 S(I,J) = S(1I,J)+C(I,J) 
115 NEXT J 
ar, 120 NEXT I 
125 N=N+1 
130 IF N>7 THEN 165 
135 FOR I= 1 TO 7 
140 FOR J= 1 TO 7 
145 B(I,J) = C(I,J) 
150 NEXT J 
155 NEXT I 
160 GOTO 75 
165 PR. "THE SUMMATION MATRIX ES :" 
170 FOR I= 1 TO 7 
| 175 FOR J= 1 TO 7 
— 180 PRINT S(I,J); 
. 185 NEXT J 
2 190 PRINT 
195 NEXT I 
200 END 
Bewal 


*Can you code this step in a more "elegant" manner? 
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Problem 2 iad 
(a) The schedule of a larger airline company is given 
below. Calculate the number of paths a passenger can 
choose in travelling from Philadelphia to Hawaii, if he ’ 
wants to make stops in any three cities before reach- 3 
ing his destination (Hawaii). = 
(b) On his return, he wishes to go from Hawaii to Washing- 
ton, again stopping in any three cities before reaching 
Washington. How many choices does the passenger now 


have? 


ANCHORAGE 


“BROT 


oe 





WINNIPEG 
MINNEAPOLIS/ST. PAUL 










YO AND / ~~ 


THE FAR EAST PORTLAND g 
PH] LADELPHIA 7 


\ PITTSBURGH = 
ia RANCISCO a ASHINGTON 
FT. LAUDERDALE 
TAMPA 
ST. PETERSBURG 
MIAMI 


LOS ANGELES 








HAWAII 


we 


Problem 3 

Communication matrices can also be used in investigating 
directed communication of messages. 

In the newly formed World Organization Opposed to Polluting 
the Sea there are eight member countries. In appointing repre- 
sentatives to WOOPS the member countries governing officers over- 
looked their appointee's language qualifications. Thus each ap- 
pointee cannot speak to all of the other representatives directly. 
It is also possible that some representatives are acting under a 
"closed-mouth" policy. This policy allows a representative to listen 
to one speaking to him but he, in turn, will not speak or respond 


to this particular representative. 


EXAMPLE: 
The Spanish representative speaks to the German representative. 


The German representative speaks to the Spanish and the Ameri- 
can representatives. 


The Italian representative speaks to the American and the French 
representatives. 


The Chinese representative speaks to the Russian representative. 


The American representative speaks to the German and the Italian 
representatives. 


The Nigerian representative speaks to the Russian representative. 


The Russian representative speaks to the Chinese, the Nigerian, 
and the French representatives. 


The French representative speaks to the Russian representative. 


(In this example the French representative is the only one oper- 


ating under the "closed-mouth" policy.) 


Finally, let's assume that a message becomes hopelessly garbled 


if it is translated more than five times. 


= £29 = 


Problem: 


The German representative wishes to argue the validity of 
throwing empty beer bottles overboard from his country's freighters. 


Can his statement be routed so that all representatives of WOOPS 


receive his message? Which countries can and which cannot get 





messages through to all representatives? 

Problem 4 - 
A strange new plague has swept over the countryside of Planet xX. 

It is due to Virus V1, which arrived via a meteorite that has land- 

ed from outer space. The inhabitants of Planet X who contacted the 

meteorite directly contract a disease called Dl. To survive they 


must receive an antidote called Al. 
NOW FOR SOME BAD NEWS: 


If an inhabitant with disease Dl contacts another inhabitant 
of X, the virus is transferred in mutated form as V2. The only 
antidote for V2 is called A2. Similarly, further transfer of 
the virus produces additional mutations with new antidotes re- 


quired as follows: 


INHABITANTS OF PLANET 


ned 
Fy “Vf 
Xe. ef”. LEAK 
." 'e s 
se 


X 


4; 






METEORITE & 
DISEASE D1 D2 D3 D4 NONE 
ANTIDOTE Al A2 A3 A4 








oe 


It is possible for an inhabitant to be exposed to several 
viruses through various contacts, in which case he must receive 
. each appropriate antidote. The amount of antidote to be given 
is proportional to the number of paths by which the virus was 


communicated, with 1 gram needed for each communication. 
AND NOW FOR SOME GOOD NEWS: 


It turns out that virus V5 (and all higher numbered viruses) 
are harmless. Thus‘only four antidotes (Al, A2, A3, and A4) are 


needed to treat any inhabitant of Planet X. 
EXAMPLE: 


Let's prescribe appropriate antidote dosages for the three 


inhabitants shown below. 


OD The arrows indicate 
[METEORITE ie communication between 
——-- inhabitants. 


mM € Oo ® ™M] ® ® @ 

a a! 0 0 2 «0 

A* =. (2) | | A? = 0 21 0 0 
! i i 

@io rt o 0° @ 0 0 1 0 

Ba Or Dev dl OY @.o0 1 0 0 

Dy 2 © iM] 2D @® © 

[mt | :_ = » 9 {M]/o0 0 2 0 

A? = : (0 0 i. 0 A*= @/o 1 0 0 

2 | 0 1 0 0 | @\o o 1 0 

@ to 0° 2 of Oie LT 8 ° 0, 
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NOTE: Only the first row, (M) should be examined in this problem, 
Since viruses can only originate with the meteorite M. 


SOLUTION : 


ds 
INHABITANT 1 1 GRAM 0 GRAMS 








RAM 
re 

INHABITANT 2 0 GRAMS 2 GRAMS 0 GRAMS 2 GRAMS 
INHABITANT 3 1 GRAM 0 GRAMS 0 GRAMS 0 GRAMS 

(CAN YOU FIND ALL THESE DISEASE COMMUNICATING PATHS?) 
NOW BACK TO ... 

PROBLEM 4: 
Prescribe drugs for the 17 inhabitants of Planet X shown on the 
cover of this module. If you don't have a computer you have our 

— 

sympathy! 
some selected answers to the above problems: 
1. You should have found a total of 18 paths from Eichelberger- 

town to Applebachsville. | 
2. (b) There are five ways to go from Hawaii to Washington with 

stops in three cities: 
Hawaii to Chicago to New York to Philadelphia to Washington 
Hawaii to Chicago to Detroit to Philadelphia to Washington 
Hawaii to Chicago to Cleveland to Pittsburgh to Washington @¢ 
Hawaii to Chicago to Cleveland to Philadelphia to Washington 
Hawaii to Seattle to New York to Philadelphia to Washington =, 
4. The dosage for Inhabitant 2 1S: 

1 gram of Al, 1 gram of A2, 2 grams of A3, 3 grams of A4. 
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